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An  Increase  in  the  critical  failure  tensile  stress,  9^,  to  8000  psi  from 
the  nominal  6500  psi  causes  a 60Z  reduction  in  target  damage,  while  a decrease 
to  5000  psi  increases  target  damage  by  over  100Z.  The  tensile  strength  of 
the  target  is  the  most  sensitive  macromechanical  physical  characteristic  in 
determining  target  damage. 

Target  damage  decreases  with  an  increase  in  Young's  Modulus  (Ej  or 
Poisson's  Ratio  (v)  and  with  a decrease  in  density 

The  approximate  relationship  of  the  macromechanical  property  (of)  to  a 
fundamental  micromechanical  property  (the  stress  intensity  factor  (Tt)  in 
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terms  of  the  numerical  calculation  cell  size  (AST  is  as  follows: 

\cWrt*  >) 
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The  following  damage  mechanisms  have  been  identified: 

o Tensile  failure  occurs  behind  the  shear  and  Rayleigh  waves, 
and  cracks  propagate  until  the  shear  wave  has  travelled 
about  1 droplet  diameter  in  the  1100  fps  impact  case. 

(~  .5  vsec  for  a 1 mm  drop). 

o The  water  drop  has  time  to  interact  with  the  surface  cracks 
it  generates  during  early  stages  of  impact. 

o The  presence  of  cracks  significantly  alters  the  stress 
fields  near  the  developing  cracks,  thereby  rendering 
elastic  stress  field  predictions  invalid  once  crack 
propagation  begins. 

o Crack  initiation  and  propagation  may  be  enhanced  by 
the  high  velocity  water  jets  which  are  generated  and 
which  interact  with  surface  roughness  and  flow  into 
surface  cracks  and  flaws. 

o The  geometry  of  the  drop  is  responsible  for  a pressure 
spike  early  in  the  impact  which  can  actuate  cracks  by 
forcing  water  into  surface  openings. 
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PREFACE 


When  an  IR  window  is  exposed  to  impacts  of  rain  droplets 
or  ice  particles  at  subsonic  velocities  (500-1100  fps) , a 
loss  in  optical  transmission  frequently  results  which  pre- 
cedes any  mass  loss  from  the  exposed  surface.  This  phenom- 
enology is  not  well  understood,  and  is  the  subject  of  this 
report.  The  authors  have  used  numerical  methods  to  examine 
the  dynamic  stresses  in  window  materials,  the  development 
of  internal  cracks  which  may  cause  transmission  loss,  and 
the  dependency  of  these  processes  on  various  basic  parameters 
characterizing  window  materials. 

The  work  was  performed  for  the  Air  Force  Materials  Labor- 
atory under  contract  number  F33615- 75-C- 5081 . The  guidance 
and  assistance  of  the  AFML  Project  Engineer,  T.  L.  Peterson, 
is  gratefully  acknowledged. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

1 . 1 Background  and  Objectives 

There  is  a need  for  substantially  improving  the  resis- 
tance of  IR  windows  to  optical  damage  resulting  from  subsonic 
impacts  of  rain  drops.  Such  damage  is  usually  referred  to 
as  subsonic  erosion,  although  it  differs  from  more  conventional 
erosion  in  that  significant  optical  damage  occurs  before  there 
is  any  mass  loss  from  the  surface.  The  apparent  reason  for 
the  optical  damage  is  internal  cracking,  which  causes  loss  of 
transmittance  as  well  as  degradation  of  optical  quality. 

Erosive  damage  of  optically  transparent  materials  by  water 
drop  impacts  has  been  experimentally  studied  (e.g.,  Reference  1) , 
and  data  have  been  obtained  regarding  the  optical  transmission 
loss  through  specific  materials  subjected  to  specific  erosive 
environments.  However,  it  is  difficult  to  generalize  from  these 
data  because  of  the  complex  nature  of  the  phenomena  involved. 

The  present  study  makes  use  of  an  efficient  numerical 
technique  for  investigating  response  of  window  materials  to 
subsonic  water  drop  impacts.  Specific  objectives  of  the  pro- 
gram were:  (1)  to  identify  the  mechanisms  responsible  for 

transmittance  loss  in  IR  windows  subjected  to  subsonic  rain 
erosion,  and  (2)  to  determine  the  sensitivity  of  these  mech- 
anisms to  basic  material  parameters. 

1 . 2 Analytical  Approach 

The  approach  taken  for  accomplishing  these  objectives 
is  based  on  numerical  code  solutions  using  the  WAVE-L  code. 
(W^'E-L  is  a 2-D  Lagrangian  finite  difference  code  based  on 
the  HEMP  code.  Reference  2.)  The  following  steps  are  involved 
in  the  analyses: 
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1.  Impact  loading  analyses  are  performed  to 
determine  the  history  of  stresses  applied 
to  a window  surface  by  water  drop  impacts. 

The  loading  stresses  are  stored  as  a func- 
tion of  position  on  the  target  surface  and 
time . 

2.  The  impact  loading  stresses  are  then  applied 
as  a surface  boundary  condition  to  investi- 
gate the  dynamic  stress  and  failure  response 
of  window  materials  to  the  impact  loading. 

The  sequential  relationship  between  these  steps  is  shown 
in  Figure  1. 

1 . 3 Summary  of  Results 

During  the  course  of  this  program  numerical  solutions 
were  obtained  which  predict  the  response  of  materials  to  water 
drop  impacts  at  675  and  1100  fps.  Detailed  analyses  of  target 
damage  sensitivity  to  elastic  and  tensile  properties,  material 
density,  impact  velocity  and  numerical  cell  size  were  performed 
for  21  cases. 

1.3.1  Impact  Loading  on  Rigid  Surface 

The  impact  dynamics  and  surface  pressure  loadings  from 
a 1 mm  water  drop  impacting  a rigid  surface  at  675  fps  and  at 
1100  fps  were  calculated  using  the  WAVE-L  code.  Figure  2 shows 
two  representative  velocity  field  plots  from  the  1100  fps  impact 
solution,  and  illustrates  the  nature  of  the  impact  process. 
Jetting,  the  development  of  radial  velocities  greater  than  the 
impact  velocity,  is  prominent  in  the  figure.  The  peak  ratio 
of  jetting  velocity  to  impact  velocity  in  the  1100  fps  numeri- 
cal solution  was  3,  which  is  consistent  with  observed  values 
(Reference  3) . 
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Figure  1.  Steps  Involved  in 
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Figure  3 shows  the  pressure-time  history  from  the  1100  fps 
sole  ion  at  the  impact  point.  The  numerical  solution  generated 
by  Huang  (Reference  4)  for  an  impact  velocity  of  980  fps  is  also 
on  this  figure.  The  solutions  are  similar  except  during  the 
first  0.1  usee,  where  the  WAVE-L  solution  shows  an  almost  in- 
stantaneous rise  to  approximately  the  water  hammer  pressure, 
in  contrast  with  the  more  gradual  rise  to  a lower  pressure 
shown  by  Huang.  However,  particular  efforts  were  made  in  the 
WAVE-L  solution  to  resolve  transitory  phenomena;  thus  we  be- 
lieve the  behavior  of  the  WAVE-L  solution  is  more  accurate 
during  the  first  0.1  ysec. 

Figure  4 presents  several  pressure-radius  profiles  from 
the  1100  fps  impact  case.  These  profiles  show,  at  several 
times,  the  spatial  variations  in  pressure  which  load  the  sur- 
face of  a window  from  a single  drop.  Note  in  the  profiles 
at  t = .12,  .24,  and  .48  ysec  that  pressures  are  greatest 
near  the  periphery  of  the  expanding  contact  surface. 

1.3.2  Window  Response  and  Sensitivity  Study  for  Water 
Drop  Impacts 

Table  1 lists  the  21  window  response  and  sensitivity 

cases  analyzed  in  this  study.  ZnSe  (p  = 5.27  gm/cm^ , 

£ ^ 

E = 10  x 10  psi,  v = .3)  was  the  baseline  window  material 
treated.  Impacts  were  simulated  both  on  purely  elastic  ZnSe 
and  on  ZnSe  which  is  allowed  to  fail  in  tension.  Two  types  of 
tensile  failure  were  modeled:  (1)  failure  whenever  a critical 

tensile  stress  (o^,)  is  exceeded  (o>o^.) , and  (2)  failure  when- 
ever a specified  time  integral,  k^,  of  the  tensile  over-stress 
is  exceeded  ( £ (o-Oj)dt>Kj,  for  o>Oj.  In  both  cases,  discrete 

tensile  cracks  can  form  and  propagate  through  the  computational 
grid.  (When  a cell  cracks,  the  WAVE-L  code  calculates  the 
orientation  and  width  of  the  crack,  and  also  properly  adjusts 
the  residual  properties  and  stresses  in  the  cell. 
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TARGET  RESPONSE  CASES  FOR  WATER  DROP  IMPACTS  (D  - . 1 cm) 
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1100  5 6500  10  .3  5.27  50  Thin  Target 


The  propagation  of  tensile  cracks  from  a 1100  fps  water 
drop  impact  is  illustrated  in  Figure  5.  This  figure  shows 
the  final  crack  pattern  from  the  Case  1 solution,  which  used 
the  critical  tensile  stress  criteria  (<T  = 0)  with  qj  = 6500 
psi.  The  tensile  cracks  propagate  at  velocities  on  the  order 
of  the  material  shear  (Cs  = .22  cm/psec)  and  Rayleigh  wave 
speeds.  In  fact,  the  cracks  form  and  propagate  behind  the 
shear  wave  front  which  develops  during  the  water  drop  impacts. 
By  .5  \isec,  the  tensile  cracks  are  completely  formed  in  this 
impact. 

The  sensitivity  study  involved  the  variation  of  material 
parameters  (kt,  aT,  E,  v,  pQ) , impact  velocity  (V),  and  compu- 
tational zone  size  (Ax).  The  crack  surface  area  (A),  as  com- 
puted by  the  numerical  solutions,  was  used  as  a measure  of  the 
target  damage.  Table  1 lists  the  target  damage  in  non-dimen- 
sional  form,  A/An,  where  An  = 8 x 10  cm  is  the  cracked  sur- 
face area  for  the  base  case  (Case  2) . 

Figure  6 summarizes  the  results  of  the  material  parameter 
sensitivity  study.  Figure  6A  shows  target  damage  dependence 
on  the  stress  integral  parameter  (tc^) . The  target  damage  is 
a relatively  strong  function  of  and  therefore  of  the  strain 

rate  characteristics  of  the  target.  An  increase  of  <T  from 

-7  * 

zero  to  5 x 10  Mbar-psec  reduces  the  target  damage  by  37%. 

_ 7 

(The  = 5 x 10  results  are  used  as  the  baseline  target 
damage  case  in  this  study.) 

Figure  6B  shows  target  damage  when  critical  tensile 
strength,  oT,  is  varied  and  = 5 x 10  Mbar-psec  is  held 
constant.  An  increase  in  to  8000  psi  from  the  nominal  Oj  = 
6500  psi  causes  a 60%  reduction  in  target  damage.  A decrease 
to  5000  psi  causes  the  target  damage  to  increase  by  over  100%. 
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Final  Crack  Pattern 


Figure  S.  Approximate  Crack  Length  versus  Time  for  1100  fps  Water 
Drop  Impact  on  ZnSe  Window  (Case  1) 
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Figure  6A,B.  Target  Damage  as  a Function  of  Stress  Integral  Parameter  (kt) 

and  Tensile  Strength  (oT) 
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The  tensile  strength  of  the  target  is  the  most  sensitive 
physical  material  characteristic  in  determining  target  damage. 

Figures  6C , 6D , and  6E  show  target  damage  caused  by 
variations  in  the  elastic  Young's  Modulus  (E)  , Poisson's  Ratio 
(v) , and  material  density  (p  ) . The  target  damage  decreases 
with  an  increase  in  E or  v>  and  decreases  with  a decrease  in 
po.  In  Section  5.3,  the  effect  on  target  damage  is  inves- 
tigated when  both  E and  pQ  are  changed.  When  E is  increased 
to  14  x 10^  psi  and,  at  the  same  time , pQ  is  decreased  to 
po  = 4 gm/cm^,  target  damage  is  decreased  by  67%. 

The  675  fps  water  drop  impact  on  the  baseline  ZnSe  target 
(K,p  = 5 x 10  Mbar-ysec,  Oj  = 6500  psi.  Ax  = .005  cm)  did  not 
generate  any  tensile  cracks  (Case  19),  since  the  peak  tensile 
stress  exceeded  Oy  only  for  a very  short  time.  The  solution 
was  repeated  with  Ky  = 0 (Case  20)  , and  one  tensile  crack 
developed.  The  final  target  conditions  were  then  used  for  a 
second  and  third  water  drop  impact  on  the  same  site.  This 
simple  multiple  impact  simulation  predicts  an  increase  in  crack 
damage. 

The  zoning  sensitivity  study  revealed  an  important 
relationship  involving  the  numerical  simulation  of  the  material 
3 tress  intensity  factor  (Kj),  which  is  a fundamental  material 
parameter  for  crack  initiation  and  propagation.  The  following 
approximate  relationship  is  found  between  Kj,  the  numerical 
zone  size  (Ax) , and  the  critical  tensile  stress  for  fracture  (Oy) : 

Kj  * (TTAX^Oj.  (1) 

This  equation  is  important  because  Kj  is  a measurable 
physical  material  parameter  which  is  more  fundamental  than  Oy 
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Figure  6C,D,E.  Target  Damage  as  a Function  of  Young's  Modulus  (E) , 
Poisson's  Ratio  (v),  and  Material  Density  (Pq) 
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from  the  micromechanics  point  of  view.  Equation  1 thus  relates 
the  macromechanical  property  (o-p)  to  a fundamental  micromechanical 
property  (KjJ  in  terms  of  the  numerical  cell  size  (ax). 

The  similarity  of  crack  patterns  and  predicted  target 
damage  from  Cases  15  and  16  (figure  7)  help  demonstrate  the  validity 
of  Equation  1.  The  Kj  for  Cases  15  and  16  are  the  same  to  within 
about  8%;  however,  o^,  decreases  by  a factor  of  = .46  while 

the  zone  size  increases  by  a factor  of  4.  The  relative  target 
damage  for  these  two  cases  is  very  similar , 3.57  (Case  15)  and 
3.36  (Case  16).  Details  of  this  material  stress  intensity 
factor  model  are  described  in  Section  5.7. 

1.3.3  Damage  Mechanisms 

The  most  important  conclusions  concerning  damage  mechanisms 
responsible  for  transmittance  loss  in  1R  windows  are  summarized 
below : 

1.  The  dynamic  tensile  stresses  which  activate  crack 
propagation  occur  behind  the  shear  and  Rayleigh 
waves  which  develop  in  the  target  during  the  water 
drop  impact.  Tensile  cracks  can  thus  propagate 
behind  the  shear  wave  until  the  tensile  stresses 
in  the  target  fail  to  satisfy  a crack  propagation 
criterion.  For  the  1100  fps  impacts  treated  in  this 
study,  the  water  drop  continues  to  load  the  target 
surface  for  over  2 psec,  but  the  crack  propagation 
lasts  only  .5  psec.  The  loading  processes  therefore 
have  time  to  interact  with  the  surface  cracks  which 
result  from  the  loading. 

2.  In  the  1100  fps  (750  mi/hr)  water  drop  impact, 

the  dynamic  tensile  stresses  initiate  and  propagate 
cracks  in  the  ZnSe  targets.  The  presence  of  the 
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Igure  7.  Crack  Patterns 


cracks  significantly  alters  the  stress  fields. 
Therefore,  elastic  stress  field  predictions  are 
not  valid  once  crack  propagation  begins.  Dynamic 
crack  propagation  problems  involve  significant 
coupling  between  the  crack  (viewed  as  a growing 
and  moving  free  surface  boundary  condition)  and 
the  non-uniformly  stressed  target  material. 

3.  The  water  drop  impacts  on  the  rigid  surface  and 
on  the  ZnSe  targets  at  675  fps  and  1100  fps  gen- 
erate high  velocity  radial  flow  of  water  along  the 
target  surface.  The  radial  velocities  calculated 
in  this  study  are  2.5  and  3 times  the  impact  velo- 
city for  the  675  fps  and  1100  fps  impacts,  re- 
spectively. The  high  velocity  jet  is  a potential 
source  of  damage  to  the  target  surface  from  (a) 
interactions  with  the  intrinsic  surface  roughness, 
and  (b)  interactions  of  high  velocity  water 
flowing  into  surface  cracks  and  flaws.  Both  of 
these  mechanisms  could  represent  important  proces- 
ses for  enhancing  crack  initiation  and  propagation 

4.  Spherical  water  drop  impacts  on  rigid  and  nearly 
rigid  surfaces  generate  an  off-axis  pressure  tran- 
sient of  about  2.7  times  the  water  hammer  pressure 
This  pressure  spike  occurs  near  the  periphery  of 
the  expanding  contact  surface,  very  early  in  the 
impact  process.  This  pressure  spike  can  actuate 
surface  cracks  and  flaws,  both  because  of  the  rela 
tively  high  stresses  which  it  induces,  and  by 
forcing  high  pressure  water  into  surface  openings. 

5.  The  target  damage  from  crack  propagation  can  be 
simulated  in  the  numerical  solutions  using  the 
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fundamental  material  stress  intensity  factor  (Kj) 
based  on  the  following  approximate  relationship. 

u 

Kj  = (ttAx)2Oj 


1 . 4 Recommendations 

The  results  of  this  study  demonstrate  that  an  efficient 
numerical  technique  for  investigating  the  response  of  window 
materials  to  subsonic  water  drop  impacts  has  been  developed. 
Also,  the  unique  capability  of  numerical  solutions  to  predict 
the  response  of  hypothetical  materials  of  possible  interest 
was  successfully  applied  in  the  form  of  a detailed  target 
damage  sensitivity  study  involving  elastic  properties,  density, 
tensile  failure  properties,  impact  velocity,  and  numerical 
cell  size. 

The  most  important  results  of  this  study  involve 
the  predictions  of  specific  crack  patterns  and  damaged  areas 
as  a function  of  tensile  failure  criteria  and  parameters.  Also, 
the  predicted  target  damage  results  based  on  the  material 
stress  intensity  factor  (Kj)  was  significant.  This  aspect  of 
the  study  should  be  continued  and  expanded  to  include  more 
detailed  examination  of  crack  propagation  theory  and  to  develop 
a model  which  connects  the  microscopic  derails  of  crack  propa- 
gation to  the  macroscopic  description  inherent  in  numerical 
codes . 

The  following  actions  are  recommended  in  order  to  take 
advantage  of  the  results  of  this  study: 

1.  Examine  and  model  the  relationship  between  micro- 
scopic crack  initiation  and  propagation  theory 
and  the  macroscopic  descriptions  which  can  be 
provided  by  numerical  codes. 
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2.  Apply  the  numerical  codes  to  other  materials  of 
interest . 

3.  Apply  the  numerical  codes  to  potential  multi- 
material window  designs  of  interest. 

4.  Use  the  numerical  codes  in  a plane  strain  mode 
to  investigate  multiple-impact  phenomenology  for 
impacts  on  separated  sites. 
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SECTION  II 


IMPACT  LOADING  ON  A RIGID  SURFACE 

The  impact  of  a water  drop  on  a surface  produces  very 
high  transient  pressures  or  stresses  in  both  the  drop  and  the 
target  material  near  the  contact  surface.  Shock  or  intense 
stress  waves  propagate  from  the  immediate  impact  interface, 
both  into  the  target  material  and  into  the  st  i 1 1 - advanc i ng  drop- 
let. Geometric  divergence  and  pressure  relief  from  free 
surfaces  result  in  rapid  reduction  of  the  peak  pressure/stress. 

The  wave  propagating  into  the  water  drop  produces  decel- 
eration as  well  as  severe  distortion  as  the  high  pressures 
interact  with  the  walls  of  the  drop.  The  drop  flattens 
against  the  target  surface,  and  the  pressure  falls  off  us  the 
region  of  impact  loading  grows.  Simultaneously,  target  mat- 
erial beneath  the  impact  is  compressed.  In  very  high-velocity 
impacts,  this  localized  compression  is  sufficient  to  produce 
significant  axial  and  radial  distortion  in  the  target,  re- 
sulting in  penetration  and/or  cratering.  In  the  relatively 
low-velocity  water  impacts  which  are  of  interest  in  this  pro- 
gram, the  target  surface  will  respond  with  only  small  distor- 
tions. This  reduces  the  complexity  of  the  impact  dynamics, 
since  loading  occurs  on  an  essentially  flat  surface. 

2 . 1 General  Procedure  for  Impact  Analysis 

In  order  to  provide  a loading  function,  P(r,t),  for  the 
analyses  of  dynamic  response  of  IR  window  materials,  two 
impact  solutions  were  calculated  for  a 1 mm  spherical  water 
drop  on  a rigid  plane  surface.  The  WAVE-L  code  was  used  for 
these  solutions. 

The  initial  water  drop  velocities  in  the  impact  problems 
are  675  fps  (460  mi/hr)  and  1100  fps  (750  mi/hr).  The  numer- 
ical fit  of  Walker  and  Sternberg  (Ref.  5)  to  the  equation  of 
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state  of  water,  as  given  in  Table  2 was  used  in  the  calcu- 
lations. Since  the  purpose  of  these  calculations  was  to  obtain 
the  pressure  on  the  plane  interface  as  a function  of  time  and 
position,  and  since  the  largest  and  potentially  most  important 
pressures  and  pressure  gradients  occur  early  in  the  impacts, 
the  following  procedure  was  adopted  for  each  solution: 

(a)  A relatively  finely  zoned  grid  (shown  in  figure  8A) 
was  used  to  determine  the  early-time  loading. 

The  vertical  axis  at  r = 0 is  the  axis  of  symmetry. 
Since  this  grid  was  designed  with  computational 
accuracy  in  mind,  the  interface  stress  function 
obtained  during  the  initial  high  stress  phase  of 
the  solution  is  well  resolved.  Integration  of  this 
grid  proceeded  until  .34  ysec  in  the  675  fps  solu- 
tions, and  until  .26  usee  in  the  1100  fps  solution. 

The  maximum  pressures  at  these  times  were  approx- 
imately 1.0  Kbar  (from  3.1  Kbar  water  hammer  pres- 
sure) and  2.6  Kbar  (from  5.1  Kbar  water  hammer 
pressure),  respectively. 

(b)  A less  finely  zoned  grid  (shown  in  Figure  8B)  was 
used  to  determine  the  pressure  loading  at  later 
times  (i.e.,  after  .34  usee  and  .26  usee  in  the 
675  and  1100  fps  impacts,  respectively.)  Since 
this  grid  was  chosen  as  a compromise  between  ease 
of  rezoning  and  computational  considerations,  the 
integration  was  not  exceptionally  difficult,  as 

it  would  have  been  using  the  grid  shown  in  Figure  8A. 

The  physical  phenomena  taking  place  during  the  impact 
of  a water  drop  have  been  examined  both  theoretically  and 
experimentally  by  many  workers  (e.g.,  References  6-9).  Generally, 
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TABLE  2. 


WALKER-STERNBERG  EQUATION  OF  STATE 
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TABLE  2-  WALKER-STERNBERG  EQUATION  OF  STATE  (Continued) 
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the  processes  have  been  described  in  two  phases.  Although 
our  numerical  solutions  make  no  distinction  between  these 
phases,  it  is  convenient  to  adopt  the  two-phase  nomenclature. 
During  Phase  1,  the  compressibility  of  the  liquid  is  the 
dominating  influence  (characterized  by  high  pressures  and 
very  little  radial  flow).  Phase  2 is  dominated  by  radial 
flow  of  the  liquid  at  velocities  which  are  generally  greater 
than  the  initial  impact  velocity. 

2 . 2 675  fps  Impact  on  Rigid  Surface 

Figures  9 and  10  show  the  velocity  field  and  pressure 
distribution  on  the  rigid  boundary  for  the  675  fps  solution 
at  .13  ysec  after  initial  impact.  At  this  time,  the  impact 
is  largely  influenced  by  the  Phase  1 compressibility  effects. 
The  drop  is  still  essentially  spherical,  but  compressed  at 
the  bottom.  There  is  very  little  radial  flow  beyond  the 
spherical  contact  surface,  and  the  peak  pressure  of  1.6  Kbar 
is  still  about  half  of  the  characteristic  water  hammer  pres- 
sure ( p0VC  =3.1  Kbar)  for  this  case.  At  this  early  time, 
the  pressure  being  applied  to  the  rigid  plane  is  highest 
toward  the  periphery  of  the  contact  area. 

Figure  11  shows  the  velocity  field  and  pressure  contours 
in  the  drop  at  .26  ysec.  Here  the  velocity  plot  shows  that 
radial  flow  is  beginning,  but  has  not  yet  caused  large  dis- 
tortions. The  pressure  is  still  high,  and  the  point  of  appli- 
cation of  the  greatest  pressure  on  the  rigid  plane  remains 
at  the  periphery  of  the  contact  area  (see  Figure  10).  It  is 
interesting  to  note  in  Figure  11  that  the  peak  pressure  region 
within  the  drop  is  not  on  the  rigid  target  surface,  but  rather 
is  slightly  above  that  surface. 
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Figure  9.  Velocity  Field  for  675  fps  Impact  of  a Water  Drop  on  a 

Rigid  Surface  at  .13  Msec 
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Figure  10.  Radial  Pressure  Profiles  on  a Rigid  Surface  for 
675  fps  Water  Drop  Impact 


Figures  12-14  show  the  progress  of  the  solutions  as 
large  radial  velocities  develop  along  the  surface  of  the 
rigid  plane  (i.e.,  the  Phase  2 processes).  The  location  of 
the  pressure  maximum  within  the  water  drop  moves  down  to  the 
rigid  interface  plane  and  the  pressure  distribution  on  the 
interface  approaches  uniformity,  with  a rapid  fall-off  at  the 
edge  of  the  contact  area  (see  Figure  10,  page  26).  The  peak 
pressure  falls  to  about  1/2  PV  (.2  Kbar  in  the  675  fps  case) 
as  the  problem  becomes  dominated  by  the  physics  of  incompres- 
sible flow.  The  loading  stresses  on  the  target  will  vanish 
as  the  rear  of  the  drop  approaches  the  target  surface  at 
about  D/V0  = 5 psec. 

The  interface  pressure  vs  radius  during  the  transition 
cycles  between  the  two  grids  used  in  the  solution  is  shown 
in  Figure  15.  Although  some  loss  of  resolution  at  the  edge 
of  the  impact  region  occurs  during  the  transition,  this  will 
have  little  effect  on  target  response  because  of  the  rela- 
tively low  stress  levels. 

2 . 3 1100  fps  Impact  Solution 

The  processes  at  work  during  the  course  of  the  1100  fps 
solution  are  fundamentally  the  same  as  those  described  for 
the  675  fps  solutions.  The  water  hammer  pressure  (pQVC)  in 
this  case  is  5.1  Kbar. 

Figure  16  shows  the  velocity  field  and  the  pressure  con- 
tours in  the  drop  at  .12  psec.  (Pressure  distributions  on 
the  rigid  surface  for  this  case  are  presented  in  Figure  4, 
page  8).  The  processes  up  to  .12  psec  are  largely  dominated 
by  compreaaible  effects.  There  is  little  grid  distortion 
(other  than  a flattening  of  the  sphere  at  the  bottom) , and 
high  pressures  remain  around  the  impact  point.  Some  radial 
flow  is  apparent  but  it  causes  no  significant  distortion  of 
the  drop. 
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Figure  16.  Velocity  and  Pressure  Fields  for  1100  fps  Water  Drop  Impact  on  Rigid  Surface  at  .12  usee 


Figure  17  shows  the  solution  at  .24  psec.  At  this  time 
Phase  2 behavior  is  apparent.  Large  radial  velocities  and 
significant  radial  distortions  of  the  drop  are  seen  as  well 
as  the  same  phenomenon  that  was  observed  in  the  675  fps 
solution;  i.e.,  the  point  of  maximum  pressure  within  the  drop 
is  slightly  above  the  rigid  plane.  Also,  the  greatest  pres- 
sures on  the  rigid  plane  continue  to  move  out  as  the  radius 
of  contact  increases  (as  shown  in  Figure  4). 

Figures  17,  18,  and  19  show  the  evolution  of  the  contin- 
uing solution.  The  water  drop  becomes  highly  distorted,  large 
radial  velocities  develop,  and  the  pressure  begins  to  decrease. 

As  in  the  previous  solution,  the  pressure  applied  to  the  rigid 
plane  approaches  uniformity  at  these  late  times. 

Figure  15,  page  32,  shows  the  interface  pressure  vs  radius 
for  the  transition  cycles  between  the  two  grids  used  in  the 
computation.  The  transition  is  reasonably  smooth  and  will 
cause  no  problems. 

2 . 4 An  Independent  Numerical  Calculation 

Huang  (Reference  4)  has  reported  a numerical  analysis  of  the 
980  fps  impact  of  a spherical  drop  on  a rigid  plane.  With 
the  exception  of  the  velocities,  the  conditions  in  his  analy- 
sis are  identical  to  the  1100  fps  problem  treated  here. 

Figure  3,  page  7 , shows  the  pressure  at  the  impact  point  as 
a function  of  time  for  the  two  solutions.  The  solutions  are 
seen  to  be  quite  similar,  except  at  very  early  times:  Huang's 

solution  shows  a slow  rise  (more  than  .1  psec)  to  a peak  pres- 
sure of  65%  of  the  water  hammer  pressure,  whereas  the  WAVE-L 
solution  shows  a virtually  instantaneous  rise  to  a peak  pres- 
sure of  115%  of  the  water  hammer  pressure.  In  the  WAVE-L 
solution,  small  numerical  cells  and  a minimum  of  artificial 
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18.  Velocity  and  Pressure  Fields  for  1100  fps  Water  Drop  Impact  on  a Rigid  Surface 


viscosity  were  used  to  allow  transitory  phenomena  to  be  re- 
solved, and  thus  we  believe  an  instantaneous  pressure  rise  to 
approximately  the  water  hammer  pressure  is  the  correct  behavior 
on  the  axis  of  cylindrical  symmetry. 

Figure  20  shows  pressure  vs  radius  plots  at  various  times 
for  the  WAVE-L  solution  and  Huang's  solution.  Again,  except 
for  the  early  time  discrepancy  noted  above,  the  two  solutions 
agree  well. 

i 

2 . 5 Experimental  Verification 

Figures  4 and  10,  pages  8 and  26,  (showing  pressure  vs 
radius  on  the  target  surface)  demonstrate  that  during  the  high 
pressure  phase  of  the  impact,  the  peak  pressures  occur  at  the 
edge  of  the  contact  area  rather  than  the  center.  Experimental 
indications  for  this  have  been  reported  by  Engle  (Reference  6)  and 
Rochester  and  Brunton  (Reference  8) . 

Jetting,  large  radial  velocities  late  in  the  impact  pro- 
cess, has  been  observed  by  ..iany  workers.  Fyall  (Reference  9)  has 
reported  ratios  of  peak  radial  velocity  to  impact  velocity 
ranging  from  2.9  to  3.7  for  2 mm  water  drops  impacting  Perspex 
at  velocities  from  904  fps  to  987  fps.  The  ratio  in  the  WAVE-L 
1100  fps  solution  was  3.0. 

2 . 6 Details  of  the  Pressurfe  Peak  at  Very  Early  Times  on  the 

Outer  Contact  Radius  Between~the  Spherical  Water  TTrop  and 

Rigid  Surface 

Both  the  675  fps  and  1100  fps  impact  cases  show  an 
interesting  and  possibly  important  characteristic  at  very  early 
stages  of  the  interactions.  As  described  in  Sections  2.2  and 
2.3,  the  loading  pressure  on  the  rigid  surface  peaks  near  the 
outer  contact  radius  during  the  early  stages  of  the  impact 
event.  This  phenomenon  was  investigated  using  a very  finely 
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zoned  grid  under  a separate  contract  with  Prototype  Development 
Associates  and  the  Office  of  Naval  Research  (Reference  10)  . A 
calculation  was  performed  to  obtain  more  detailed  information 
on  the  target  pressure  load  at  very  early  times  for  a very 
finely  resolved  water  drop  impact  at  1100  .'ps.  In  this  study, 
a peak  pressure  of  2.65  times  the  water  hammer  pressure  is 
predicted.  Figure  21  illustrates  the  time  resolved  pressure 
spike  at  several  radii. 

The  general  nature  of  the  off-axis  pressure  peak  can  be 
understood  in  terms  of  the  parameters  in  the  following  sketch: 


V = impact  velocity 
r * contact  radius 

R-Vt  = distance  from  water  drop  center 
to  rigid  surface  at  time  t 

The  following  relationship  is  apparent: 

r2  + (R-Vt)2  = R2  (2) 

or  in  non-dimensional  form 

(J)2  + (1  - x)2  " 1 (3) 
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Equation  3 represents  the  relationship  between  the  contact  radius 

(r)  and  time  after  impact  (t)  . Taking  the  time  derivative  of 

Equation  3 and  solving  for  r (where  r is  the  speed  of  the 

V' 

contact  radius),  we  find 


r 


V 


1 - Vt 

l 1C 

r/R 


(4) 


Thus , Equation  4 combined  with  Equation  3 relates  r to  time  or  radius. 

V 

Figure  22  is  a graph  showing  the  relationship  between  r/R, 

Vt  r 

, and  — . Note  that  because  of  the  spherical  geometry  of 

the  drop,  r is  infinite  at  t * 0.  The  results  of  this  study 
and  Reference  10  indicate  that  the  pressure  peak  at  the  contact 
radius  (r)  increases  with  time  until  the  speed  of  the  contact 
radius  (r)  equals  the  local  sound  speed  or  shock  speed  (C)  of 
the  fluid  drop. 

When  the  contact  radius  velocity  falls  below  the  local 
sound  speed  in  the  drop,  i.e.,  r<C,  the  peak  pressure  on  the 
contact  radius  begins  to  decay  from  interactions  with  the  near- 
by free  surfaces.  As  long  as  the  contact  radius  is  moving 
faster  than  any  wave  velocity,  the  drop  free  surfaces  provide 
no  relief  to  the  pressure  build-up.  Table  3 and  Figure  22 

Vt  f 

show  the  critical  values  of  r/R,  , and  y-  when  r = C = 5000 

fps  for  the  two  impact  cases  treated  in  this  study,  i.e., 

V = 675  fps  and  1100  fps. 
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TABLE  3. 

CRITICAL  WATER  DROP  PARAMETERS  WHEN  r = 
5000  FPS  FOR  V = 675  AND  1100  FPS 

C = 

V = C = 675  fps 
(0.020514  cm/psec) 

V = 1100  fps 
(0.033528  cm/psec) 

r/V 

7.41 

4.55 

r/R 

0.13 

0.22 

Vt/R 

0.009 

(t=0.022  psec  for 
R=0 . 05  cm  ) 

0.023 

(t=0.034  psec 
R=0.05  cm) 

for 

Thus,  the  peak  pressure  spike  occurs  early  in  the  drop- 
target  interaction  process  and  at  radii  which  are  only  13% 
and  22%  of  the  drop  radius  for  675  and  1100  fps  water  drop 
impacts,  respectively.  The  predicted  time  and  location  of 
the  pressure  spike  agrees  very  well  with  the  results  shown 
in  Figure  21. 

2 . 7 Comparison  of  Impacts  of  a Water  Drop  on  a Rigid  Target 
and  on  a Znse  Window 

When  the  pressure  applied  to  the  window  by  the  impacting 

water  drop  begins  to  decay,  relief  waves  propagate  into  the 

window  material  and  the  window  "rebounds".  As  a validation 

* 

of  the  assumption  that  water  impact  loading  on  a rigid  surface 
is  essentially  the  same  as  water  impact  loading  on  an  IR  window 
(in  the  velocity  range  of  interest)  a solution  of  a water  drop 
impacting  directly  on  a window  surface  was  obtained. 

Figure  23  shows  the  initial  conditions  for  this  solution. 
The  grid  design  in  the  water  drop  is  identical  to  that  used 
at  early  times  in  the  rigid  surface  solution  (i.e.,  Figure  8A) . 
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The  initial  velocity  of  the  drop  is  1100  fps.  The  window  is 
.32  cm  thick  ZnSe,  with  64  cells  across  the  thickness.  The 
two  grid  lines  which  comprise  the  interface  between  the  water 
drop  and  the  window  are  frictionless  slide  lines,  which  means 
they  can  freely  slip  along  one  another.  In  the  radial  direction 
the  grid  extends  to  .6  cm  (the  entire  grid  is  not  shown  in 
Figure  23) . 

Figure  24  compares  the  velocity  field  from  this  solution 
at  .26  ysec  with  that  of  the  1100  fps  rigid  surface  impact 
described  previously.  They  are  almost  identical.  Figure  25 
compares  the  interface  pressure  functions  for  the  two  solu- 
tions at  .26  usee,  and  except  for  numerical  oscillations,  there 
is  no  substantial  difference. 

By  .26  usee  the  peak  pressures  in  the  1100  fps  impact 
case  have  decayed  from  almost  6 Kbar  to  less  than  2.5  Kbar, 
and  the  window  material  at  the  impact  point  has  begun  to 
rebound  after  negligible  displacement.  If  the  deformation  of 
the  window  were  going  to  affect  the  interface  loading,  this 
would  have  become  evident  by  this  time.  Consequently,  we 
conclude  that  the  rigid  surface  concept  is  valid  in  the  1100 
fps  velocity  regime,  and  that  the  impact  pressure  function 
determined  using  a rigid  surface  can  be  used  to  load  elastic 
windows  without  alteration. 
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Figure  23.  Initial  Grid  Configuration  for  1100  fps  Impact  of  Water  Drop 
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Figure  25.  Pressure-Radius  Profiles  on  the  Impact  Surface  for 
1100  fps  Impact  of  a Water  Drop  on  a Rigid  Surface 
and  a ZnSe  Window 
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SECTION  III 


MATERIAL  PROPERTIES 

The  baseline  material  for  this  study  is  ZnSe  in  a form 
suitable  for  the  construction  of  optical  components  from  the 
Eastman  Kodak  Company,  under  the  name  IRTRAN  4,  and  from  the 
Raytheon  Corporation.  IRTRAN  4 is  made  by  subjecting  a ZnSe 
powder  to  heat  and  pressure  while  bulk  diffusion  and  plastic 
flow  take  place  and  a glass-like  form  of  the  material  is 
produced.  The  Raytheon  product  is  produced  by  a chemical 
vapor  deposition  process.  References  11-18  were  reviewed  to 
determine  the  properties  of  the  various  types  of  available 
ZnSe.  In  addition,  pertinent  information  was  obtained  by 
private  communications  with  J.  Pappis  of  the  Raytheon  Corp- 
oration; Roy  Rice  of  the  Naval  Research  Lab;  John  Wurst  of  the 
University  of  Dayton;  and  A.  G.  Evans  of  the  Rockwell  Science 
Center.  The  results  of  these  inquiries  are  summarized  in 
Table  4. 

ZnSe  is  modeled  in  the  computer  solutions  as  a linear 
elastic  material  with  Young's  Modulus  and  Poisson's  Ratio  as 
given  in  Table  4.  (These  imply  a bulk  modulus  K = 574  Kbar.) 
Plastic  behavior  has  not  been  observed  in  ZnSe,  though  experi- 
mental data  are  available  only  for  failure  in  tension.  We 
have  assumed  that  at  the  stress  levels  of  interest  in  the 
present  effort  compressive  yielding  does  not  occur.  Fracture 
is  allowed  to  occur,  with  the  width  and  trajectory  of  any  crack 
that  forms  being  determined  by  the  local  stress  field.  The 
stress  field  itself  will  be  appropriately  affected  by  the 
existence  of  the  fracture. 

As  a port  of  this  effort,  the  properties  of  calcium  alumino- 
silicate glass  were  compared  with  ZnSe  to  determine  if  it 
would  respond  similarly  to  water  drop  impact.  A positive 
finding  would  allow  the  use  of  calcium  a lumino - s i 1 icate  glass 
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PROPERTIES  OF  ZnSe  and  CALCIUM  ALUMINO- SILICATE  GLASS 
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Stress  Intensity  600-1300  psi  \TEn 

Factory,  K_ 


in  experimental  comparisons,  since  glass  is  more  economical 
and  has  superior  optical  characteristics.  The  particular  glass 
chosen  for  the  comparison  was  Coming's  CORTRAN  9753  , whose 
properties  are  listed  in  Table  4. 

Comparing  the  values  in  Table  4,  we  note  the  following: 

o The  density  of  ZnSe  is  almost  twice  that  of  CORTRAN 
9753  glass.  This  is  a serious  dissimilarity,  since 
it  means  that  any  inertial  effects  will  be  substan- 
tially different. 

o Both  the  Knoop  hardness  and  the  modulus  of  rupture 
of  the  9753  glass  are  greater  than  those  of  ZnSe , 
indicating  a higher  fracture  resistance  in  the  glass. 

o The  sound  speed  in  9753  glass  is  about  60%  greater 
than  that  of  ZnSe.  This  would  lead  to  differences 
in  the  timing  of  the  various  stress  wave  interactions 
(though  these  might  be  remedied  by  scaling  the  window 
thickness  and/or  the  drop  size  accordingly) . 

In  addition,  Schmitt  (Reference  1)  has  obtained  experimental 
results  in  which  both  ZnSe  and  CORTRAN  9753  glass  were  sub- 
jected to  erosion  by  water  drop  impact.  Photographs  of  the 
eroded  surfaces  indicate  that  significantly  different  processes 
are  governing  the  erosion  of  the  two  materials. 

The  property  differences  between  ZnSe  and  CORTRAN  9753, 
together  with  Adler's  data,  suggest  that  ZnSe  should  be  used 
instead  of  9753  glass  in  any  experimental  comparisons  despite 
the  associated  costs. 
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SECTION  IV 


PHENOMENOLOGY  OF  TARGET  RESPONSE  TO  WATER  DROP  IMPACTS 

All  of  the  target  response  problems  were  run  using  the 
WAVE-L  computer  code,  and  all  used  the  impact  loading  functions 
described  in  Sections  2.2  and  2.3.  The  water  drop  is  repre- 
sented as  a pressure-boundary  condition  on  the  ZnSe  plate  and 
does  not  explicitly  appear  in  the  computational  grid.  The 
numerical  solutions  therefore  become  easier  and  less  expensive. 

ZnSe  in  the  solutions  was  modeled  as  an  elastic  material 
(implying  Young's  Modulus  E = 10  x 10^  psi  and  Poisson's  Ratio 
v = .3,  with  bulk  modulus  K = 574  Kbar,  shear  modulus  G = 265 
Kbar)  but  the  ZnSe  was  allowed  to  fail  by  undergoing  tensile 
fracture  when  the  maximum  principal  stress  or  stress  history 
in  any  cell  exceeded  some  failure  criterion. 

The  WAVE-L  code  has  the  capacity  to  allow  the  formation 
of  tensile  cracks,  to  calculate  the  width  of  the  cracks,  to 
make  appropriate  changes  in  residual  properties  of  adjacent 
material,  and  to  appropriately  adjust  stresses  adjacent  to 
the  cracks.  After  a cell  has  cracked,  no  tensile  stress 
across  the  crack  is  permitted,  and  no  shear  stresses  are  per- 
mitted on  the  crack  if  the  crack  is  open  (the  width  of  the 
crack  is  continuously  monitored).  These  stress  adjustments 
are  important  because  a crack  frequently  alters  the  local 
stress  field  in  such  a way  as  to  enhance  its  own  growth;  thus, 
any  realistic  method  of  predicting  crack  growth  must  consider 
the  altered  stress  field.  Each  computational  cell  is  allowed 
to  develop  up  to  two  cracks.  If  a third  crack  forms,  the  cell 
is  considered  to  be  completely  shattered , and  no  tension  at 
all  is  thereafter  allowed.  The  general  treatment  of  tensile 
failure  and  the  plotting  conventions  used  to  indicate  the 
cracked  cells  in  the  computational  grid  plots  are  illustrated 
in  Figure  26. 
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IXJ 

F igure 


Tensile  stresses  (a)  exceed  the 
'critical  tensile  failure  criterion 
in  the  indicated  principal  stress 
direction • 

Plane  which  experiences  tensile 
failure . 


Crack  opens  and  relieves  the 
tensile  stress  such  that  no  forces 
or  tractions  exist  on  the  crack 
surfaces;  i.e.,  zero  normal  and 
shear  stress  on  the  crack  surface. 


In  the  directions  perpendicular 
to  the  crack  surface,  the  material 
can  maintain  tensile  stresses 
until  the  critical  tensile  failure 
criterion  is  exceeded  and  a new 
crack  is  formed. 


Plotting  Convention: 

Single  crack  in  r-z  plane.  The  line  is  plotted 
in  the  direction  of  the  crack  (normal  to  the 
principal  stress  causing  the  crack). 


Hoop  crack;  i.e.,  °ee  > °c  caused  the  crack. 


Hoop  and  r-z  plane  crack. 


Two  r-z  plane  cracks;  only  hoop  tension  is 
permitted . 


Completely  fractured  - no  tension  permitted. 

26.  Treatment  of  Tensile  Cracks  in  WAVE-L  Code 
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Case  1,  an  1100  fps  impact  using  the  simple  critical  stress 
failure  criterion  (see  Section  5.1  ) with  Oj  = 6500  psi, 
is  used  to  illustrate  the  phenomena  involved  in  the  target 
response  to  subsonic  water  drop  impact.  All  the  target  plots 
shown  in  the  report  use  the  target  orientation  shown  in  the 


The  water  drop  is  generally  not  shown,  but  all  impacts  are 
from  the  left  as  illustrated. 

Figure  27  shows  the  evolution  of  the  cracks  in  Case  1 
as  the  solution  progresses.  The  radial  cracks  are  initiated 
on  the  surface  at  a radius  of  about  .025  cm  (one-half  the 
radius  of  the  impacting  drop)  and  extend  outward  at  an  angle 
of  about  45°  to  the  surface  of  the  window.  This  is  consis- 
tent with  Adler's  (Reference  19),  experimental  observations  for  water 
drop  impacts  on  ZnSe,  which  show  a distribution  of  concentric 
ring  fractures. 

The  location  of  the  shear  wave  fronts  (Cg  = .224  cm/psec) 
in  Figure  27  suggests  that  the  tensile  cracks  form  and  propa- 
gate behind  the  shear  wave  front.  Figure  5,  page  10,  identifies 
two  discrete  cracks  in  the  final  crack  pattern  for  Case  1 and 
shows  crack  length  vs  time  for  these  cracks.  These  results 
demonstrate  that  the  cracks  propagate  at  a major  fraction  of 
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Figure  27. 


Crack  configuration  at  Various  times  from  Numerical  Solution  of  1100  fps 
Water  Impact  on  ZnSe  Window 


the  shear  wave  velocity.  The  formation  of  cracks  during  sub- 
sonic water  drop  impact  for  the  cases  treated  is  essentially 
over  within  the  time  required  for  the  shear  wave  to  propagate 
about  two  drop  radii. 

Further  evidence  that  the  tensile  stresses  causing  failure 
in  the  target  develop  behind  the  shear  wave  front  is  illustrated 
in  Figures  28  and  29,  which  show  contours  of  instantaneous 
peak  tensile  stress  for  Case  1 at  .2,  .3,  .4,  and  .5  ysec.  The 
shear  wave  front  is  also  shown  on  these  figures.  The  target 
volume  experiencing  some  tension  is  approximately  bounded  by 
the  shear  wave  front.  Note  that  the  Rayleigh  wave  front  is 
at  nearly  the  same  location,  since  the  Rayleigh  wave  speed 
is  = .208  cm/  ysec  or  93%  of  the  shear  wave  speed.  The 
dilatational  wave  front  is  much  farther  out  since  this  wave 
speed  is  = .42  cm/ysec  or  1.9  times  faster  than  the  shear 
wave  speed. 

A comparison  of  the  location  of  the  tensile  stress  con- 
tours in  Figures  28  and  29  with  the  crack  patterns  in  Figure  27 
shows  that  a macroscopic  tensile  stress  concentration  of  approx- 
imately the  tensile  strength  of  the  target  (o^  = 6500  psi  = 

.45  Kbar)  forms  at  the  end  of  the  running  crack.  This  compari- 
son is  made  in  Figure  30  for  t = .2  ysec.  The  crack  propa- 
gation characteristics  are  determined  by  the  tensile  stress 
concentration  near  crack  tips,  the  tensile  stress  interactions 
with  the  impact  stress  loads,  and  the  crack-free  surface 
boundary  condition. 

Figure  31  shows  contours  of  the  maximum  tensile  stresses 
experienced  in  a purely  elastic  target  (Case  17).  The  Case  1 
target,  which  fails  at  a tensile  stress  of  6500  psi  = .45  Kbar, 
shows  significant  (>  .4  Kbar)  tensions  below  the  surface.  In 
the  elastic  case,  on  the  other  hand,  high  tensile  stresses  are 
present  only  near  the  target  surface. 
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Figure  28.  Contours  of  Peak  Tensile  Principal  Stress  for  1100  fps 
Water  Drop  Impact  on  ZnSe  Window  at  .20  and  .30  psec 

7 


5 


s 

CJ 


(0 

3 

•H 

CO 

& 


0 .05 

Depth  (cm) 


Computational  Grid  Plot 
Showing  Cracked  Cells 


t = .2  psec 


Contours  of  Most- 
Tensile  Principal  Stress 


Figure  30.  Computational  Grid  and  Peak  Tensile  Principal  Stress  Contours 
for  1100  fps  Water  Drop  Impact  on  ZnSe  Window  at  .2  psec 
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SECTION  V 


SENSITIVITY  STUDIES  ON  VARIOUS  PARAMETERS  INVOLVED 
IN  WATER  DROP  IMPACTS  ON  WINDOWS 

The  1100  fps  water  drop  impact  was  used  in  all  the  sen- 
sitivity studies  except  for  the  cases  where  variation  in  impact 
velocity  was  being  examined  (Section  5.5). 

5 . 1 Tensile  Failure  Criteria  and  Parameters 

Two  tensile  failure  criteria  were  used  in  the  targets 
with  the  crack  initiation  and  propagation  model  described  in 
Section  4.  The  two  models  are  described  below: 


o Critical  Tensile  Stress  Criterion 

A tensile  crack  is  formed  whenever  the 
tensile  stress  in  the  material  exceeds  a pre- 
specified critical  stress  value,  Oj.  Thus 
failure  occurs  in  a computational  cell  if 

a>Oj . 

o Stress  Integral  Failure  Criterion 


A tensile  crack  is  formed  whenever  the 
following  integral  criterion  is  satisfied: 


l 


(O  - Oj)dt>Ky 


where 


(o  - 

O J ) 

= (o  - oT) 

if  o i.  Oj  and 

(o  - 

O-p  ) 

= 0 if  o<Oy 

0 

most 

tensile  pr 

incipal  stress 

0»r 

cr  it 

ical  stress 

value;  normally  the 

1 

stat 

ic  tensile 

strength 

Kj  * 

stre 

ss  integral 

failure  parameter 

t = time  since  stress  loading  began 
(Note  that  when  = 0,  the  stress  integral  criterion 
becomes  exactly  equivalent  to  the  critical  stress  criterion 


) 
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Figure  5,  page  in,  and  Figure  27,  page  55,  show  the 
dynamic  target  response  in  terms  of  crack  patterns  using  the 
simple  critical  stress  criterion  (<j  = 0)  with  o,p  = 6500  psi 
for  the  1100  fps  impact.  Directly  under  and  near  the  impact 
point,  cracks  form  within  the  first  few  cycles  of  the  numeri- 
cal calculation.  These  quickly  close  as  tne  surface  impact 
pressure  develops.  This  is  inconsistent  with  experimental 
observations  and  with  physical  intuition.  Tensile  stresses 
may  initially  occur  near  the  impact  point,  but  their  dura- 
tion is  so  short  that  there  probably  is  not  sufficient  time 
for  a crack  to  form.  In  order  to  investigate  the  importance 
of  tensile  stress  duration  on  crack  initiation  and  propaga- 
tion, the  stress  integral  failure  parameter,  , was  varied 
This  model  has  been  examined  for  metals  in  Reference  20. 


Figure  32  shows  the  final  crack  patterns  using  various 
values  of  Kj  for  Oy  = 6500  psi  and  for  the  1100  fps  impact. 
A value  of  Kj  = 5 x 10'?  Mbar-psec  is  sufficient  to  remove 
the  very  short  duration  cracks  near  the  impact  point.  A 
of  5 x 10'^  Mbar-psec  significantly  reduces  the  extent  of 
crack  propagation.  Figure  32  illustrates  the  qualitative 
effect  of  Kj  variations  on  target  damage. 
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The  quantitative  measure  of  target  damage  adopted  in 
this  study  is  based  on  cracked  surface  area  available  for  opti- 
cal degradation.  hhen  a computational  cell  fails  due  to  ten- 
sion in  the  target,  a cylindr i cally  symmetric  surface  is 
formed,  as  illustrated  in  the  following  sketch: 


Computational  Grid 
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Physical  Target 
(Top  view) 


RROIUS  (CM)  RROIUS  (CM) 


\ 


X (CM) 


Case  3 Case  4 

- 6 —6  . 
k = 2.0  x 10  Mbar-usec  k = 5.0  x 10  Mbar-psec 


Figure  32.  Crack  Patterns  a6  a Function  of  Stress  Integral  Parameters 
for  1100  f ps  Impact  of  . 1 cm  Water  Drop  on  ZnSc  Window 
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Optical  degradation  is  presumably  proportional  to  the 
projected  area  ot  cracked  surface,  as  shown  in  the  top  view 
in  the  sketch.  This  area,  summed  over  all  the  cracked  cells 
in  the  final  computed  crack  pattern,  is  therefore  a quanti- 
tative measure  of  target  damage.  In  order  to  compare  damage 
in  different  impacts,  it  is  convenient  to  normalize  the  calcu- 
lated projected  crack  areas  by  the  crack  area  -in  the  baseline 

. 7 

case.  The  baseline  impact  is  Case  2,  which  used  = 5 x 10 
Mbar-ysec  and  o^,  = 6500  psi  with  the  stress  integral  failure 
criterion.  For  this  case,  the  two  cracks  shown  in  Figure  32 
generate  a target  damage  surface  area  (projected  on  the  tar- 
get  surface)  of  An  = 8.07  x 10  cm  , which  is  3%  greater  than 
the  projected  area  of  the  water  drop. 

Once  this  quantitative  measure  is  established,  the 
variation  of  target  damage  for  any  sensitivity  calculation 
can  always  be  described  as  the  ratio  of  the  damaged  area  in 
the  sensitivity  case  to  the  damaged  area  in  the  baseline  case, 
A 

i.e.,  t—  . Figure  6A,  page  11,  and  Table  I,  page  6 , illus- 
An 

trate  the  variation  of  target  damage  with  based  on  the  re- 
sults shown  in  Figure  32. 

Holding  constant  the  stress  integral  failure  critei ion 
parameter  of  = 5 x 10  Mbar-ysec,  the  influence  of  the 
static  tensile  strength  parameter  (o^,)  was  then  investigated. 
Figure  33  shows  the  influence  of  static  tensile  strength  on 
predicted  final  crack  patterns  for  Oj  = 5000,  6500,  and  8000 
psi.  This  +.1500  psi  variation  of  Oj  (from  the  baseline  6500 
psi  tensile  strength)  caused  significant  changes  in  the  ex- 
tent of  crack  propagation.  For  Oj  = 5000  psi,  the  degree  of 
target  cracking  increased  significantly,  A/An  = 2.26.  The 
O-j,  = 8000  psi  target,  on  the  other  hand,  showed  only  a limited 
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Figure  33.  Crack  Patterns  as  a Function  of  Tensile  Strength,  aT,  for  1100  fps  Impact  of 

.1  cm  Water  Drop  on  ZnSe  Window 


amount  of  crack  propagation,  A/An  = .4.  Thus,  even  though 
the  purely  elastic  target  solution  (discussed  below)  showed 
tensile  stresses  up  to  about  1.6  Kbar  (23,000  psi),  a modest 
increase  of  strength  from  6500  psi  to  8000  psi  is  very  effec- 
tive in  limiting  the  crack  growth.  Figure  6B,  page  11  , illus- 
trates the  variation  in  target  damage  due  to  variations  in 
static  tensile  strength. 

Since  increasing  or  is  very  effective  in  reducing 
target  damage,  how  high  must  and  become  to  prevent  ten- 
sile cracks  from  forming?  A finely  zoned  purely  elastic  solu- 
tion was  performed  (Case  18)  to  help  answer  this  question. 

Figure  34  shows  the  peak  values  of  tensile  stress  as  a 
function  of  radius  in  the  target.  A peak  tensile  stress  of 
almost  1.6  Kbar  (23,000  psi)  is  experienced  near  the  surface 
at  about  one-half  the  water  drop  radius.  Note  that  the 
radial  gradient  of  peak  tensile  stress  is  very  sharp  as  the 
maximum  1.6  Kbar  level  is  approached. 

Figure  34  also  shows  the  peak  values  of  stress  integral 

| (o  - oT)dt)  with  oT  = 6500  psi  as  a function  of  radius  in 

the  target.  The  stress  integral  is  nearly  zero  for  radii  up 
to  about  .015  cm  ( 3R) . Then  the  stress  integral  rapidly 
peaks  to  a maximum  value  of  about  3.5  x 10'^  Mbar-ysec. 

5 . 2 Elastic  Property  Variations 

The  nominal  values  of  the  Young's  Modulus  (E)  and  Poisson's 
Ratio  (v)  for  ZnSe  are 

E = 10  x 10^  psi  x .689  Mbar 
v * .3 
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Figure  35  shows  the  crack  patterns  caused  by  varying 
Young's  Modulus  by  ±40%  from  the  nominal  value.  Crack  form- 
ation and  propagation  decreases  with  increasing  Young’s  Modulus. 
Figure  6C , page  13,  shows  target  damage  as  a function  of  vari- 
ation in  Young's  Modulus. 

Figure  36  shows  the  crack  pattern  caused  by  varying 
Poisson's  Ratio  (v  = .2,  .3,  and  .4).  Crack  damage  decreases 
as  v increases . The  nature  of  the  final  crack  pattern  is  af- 
fected by  Poisson's  Ratio.  For  v = -2,  the  cracks  begin 
closer  to  the  axis  and  propagate  much  deeper  than  for  v = -4. 
Figure  6D,  page  13,  shows  the  variation  in  target  damage  as 
a function  of  Poisson's  Ratio. 

5 . 3 Variation  of  Target  Density 

The  baseline  ZnSe  initial  density  of  pQ  = 5.27  gm/cm^ 
was  varied  in  two  calculations.  Figure  37  summarizes  the 

3 

results  of  the  calculations  for  pQ  = 4,  5.27  and  6.5  gm/cm 
in  terms  of  the  final  crack  patterns  from  the  1100  fps  water 
impacts.  Crack  damage  decreases  with  a decrease  in  initial 
density.  Figure  6E , page  13,  shows  the  variation  in  target 
damage  as  a function  of  initial  density. 

5 . 4 Variation  of  Density  and  Young’s  Modulus 

The  fact  that  increasing  Young's  Modulus  (E)  or  decreasing 
initial  density  (PQ)  limits  the  growth  of  cracks  leads  to  the 
question  of  how  effective  two  simultaneous  parameter  varia- 
tions could  be  in  reducing  crack  propagation.  Figure  38  shows 
the  results  of  calculations  in  which  PQ  and  E were  changed 
from  their  baseline  values,  separately  and  simultaneously. 

The  reduction  in  target  crack  damage  when  both  E and  PQ  are 
changed  (A/An  = .33)  is  much  greater  than  when  either  E (A/An=  .64) 
or  (A/An  = .77)  is  independently  changed. 
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Figure  38.  Crack  Patterns  as  a Function  of  Young's  Modulus  and 
Density  for  1100  fps  Impact  of  . 1 cm  Water  Drop  on 
ZnSe  Window 


5.5 


675  fps  Impact  Velocity  Case  and  Repeated  Impact  on 
the  Same  Site 

The  baseline  velocity  in  the  sensitivity  study  was  1100 

fps.  The  675  fps  impact  loading  pressure  function  was  used 

to  examine  the  influence  of  impact  velocity  on  target  damage. 

. 7 

Using  the  baseline  target  parameters,  i.e.,  = 5 x 10  Mbar- 

ysec,  o^,  = 6500  psi,  and  Ax  = .005  cm,  the  ZnSe  target  did  not 

crack  in  a 675  fps  impact  (Case  19).  Figure  39  shows  the  peak 

tensile  stress  in  the  target  as  a function  of  radius  for  this 

case.  The  maximum  tensile  stress  in  the  solution  is  31 

above  the  static  tensile  strength  (aT)  of  the  material,  but 

the  duration  of  this  overstress  is  far  too  short  to  cause 

. 7 

failure  when  k^,  = 5 x 10  Mbar-ysec. 

The  numerical  solution  was  repeated  using  the  simple 
critical  tensile  stress  fracture  criteria  (k^  = 0)  with  a^,  = 

6500  psi,  (Case  20).  In  this  case,  only  one  cell  cracked. 

Thus,  the  numerical  solution  simulated  a crack  initiation  with 
no  crack  propagation.  This  situation  is  somewhat  analogous 
to  a single  water  drop  impact  which  causes  only  marginal  damge. 
Target  damage  will  develop  after  multiple  impacts. 

Generally,  multiple  impacts  will  not  occur  on  exactly 
the  same  site.  However,  in  an  axisymmetric  code,  multiple 
impact  loadings  can  only  be  simulated  for  impacts  on  the  same 
site.  The  multiple  impact  case  on  the  same  site  is  nevertheless 
of  interest  because  of  the  interactions  between  stress  waves 
and  pre-existing  cracks.  To  investigate  the  nature  of  crack 
growth  from  stress  interactions  during  multiple  impacts  on 
the  same  site,  a second  water  drop  impact  and  then  a third 
were  simulated  on  the  target.  The  second  and  third  drop  were 
imulated  using  the  675  fps  loading  function  applied  to  the 
i il  damaged  target  configuration  from  the  previous  impact. 
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Figure  39.  Peak  Tensile  Stress  as  a Function  of  Radius 
for  675  fps  Impact  Case 
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Figure  40  shows  the  crack  pattern  after  each 
impact.  The  effect  of  the  multiple  impacts  is 
tensile  failure.  Also,  the  crack  pattern  from 
mulitple  impact  loading  appears  shallower  than 
single  impact  cases. 


water  drop 
to  spread 
the  675  fps 
in  the  1100 


fps 


5 . 6 Target  Thickness 

The  baseline  target  thickness  in  this  study  was  1/8”  = 
.3175  cm.  At  this  thickness,  the  dilatational  wave  takes 
about  1.5  psec  to  reflect  off  the  target  rear  surface  and 
interact  with  the  tensile  cracks  which  develop  near  the 
impact  site.  Since  the  1100  fps  baseline  calculation  was 
numerically  integrated  to  over  4 usee,  while  the  crack  pattern 
did  not  change  after  .5  usee,  the  1/8"  target  acted  as  if  it 
were  semi  - inf inite . 


A much  thinner  target  was  simulated  (Case  21)  to  in- 
vestigate the  potential  interaction  between  target  free  sur- 
faces and  crack  propagation.  In  Case  21,  a 1 mm  thick  target 
with  baseline  ZnSe  characteristics  was  impacted  at  1100  fps 
with  the  1 mm  water  drop.  Figure  41  shows  the  final  crack 
pattern  for  this  thin  target  case.  Crack  segments  which  formed 
in  the  .1  cm  target  but  not  i ^ ^e  .3175  cm  target  are  shown 
circled. 

The  thin  target  permittf'.'^  suff  ic  ient  interactions  with 
the  rear  surface  to  somewhat  enhance  crack  propagation.  Also, 
the  rear  of  the  target  experienced  failure  from  hoop  and 
radial  tensile  stresses.  Normally,  we  would  expect  a spall- 
type  failure  in  thin  plates,  i.e.,  from  axial  tensions.  In 
this  case,  the  fact  that  the  loading  function  has  a peak  com- 
pressive load  which  occurs  off  axis  causes  relatively  large 
radial  motions  and  stresses  which  are  larger  than  the  axial 
motions  and  stresses  at  the  target  rear  surface. 
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5.7 


Variations  in  Computational  Zone  Size  and  effective 
Va r iat ion~Tn  Target  toughness 

Figure  42  shows  effects  of  zoning  variations  in  the  final 
crack  patterns  for  three  1100  fps  water  drop  impacts  on  the 
baseline  ZnSe.  The  simple  critical  stress  failure  criterion 
was  used,  with  = 6500  psi.  The  computational  cell  or  zone 
size  (Ax)  in  the  three  cases  is  .01  , .1)05  , .0025  cm  or  tv  = 

10,  20,  and  40  (where  DQ  - .1  cm,  the  water  drop  diameter). 
Figure  42  shows  that  the  crack  propagation  characteristics 
depend  on  the  zone  size  used  in  the  calculation. 


The  stong  dependence  of  crack  propagation  on  zone  size 
is  initially  disturbing  from  the  point  of  view  of  the  ability 
of  numerical  codes  to  predict  crack  propagation.  However, 
once  the  nature  of  the  numerical  crack  propagation  model  is 
examined,  we  find  that  the  numerical  zone  size  is  related  to 
an  important  physical  material  parameter , namely  the  stress 
intensity  factor. 


The  stress  intensity  factor,  Kj,  is  a fundamental  mater- 
ial parameter  involved  in  predictions  of  crack  initiation  and 
propagation.  The  Griffith  Crack  Theory  for  linearly  elastic 
plates  predicts  the  following  relationship  for  crack  initi- 
ation in  plane  stress  or  plane  strain  configuration: 


\/~nC 


(plane  stress)  , 


°T 


I 


t: 


(i  V)c 


(plane  strain)  (5) 


where 

Kj  = stress  intensity  factor  or  crack  toughness 

2C  = length  of  pre-existing  crack  nucleus 

oj  = critical  far - f icl d tens i 1 e stress  which 
will  lead  to  unstable  crack  growth 
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Figure  42.  Crack  Pattern  for  1100  fps  Water  Drop  Impact  on  ZnSe  Window,  Zone  Size.  D/Ax 


Equation  5 predicts  the  critical  tensile  stress  for  unstable 
crack  growth  for  a uniformly  applied  tensile  stress  Oj  with 
a pre-existing  crack  of  length  2C . In  the  numerical  solution, 
crack  initiation  was  not  considered  to  be  a function  of  2C . 
Furthermore,  since  the  intrinsic  distribution  of  pre-existing 
crack  lengths  in  the  material  is  independent  of  numerical  zone 
size.  Equation  5 cannot  be  used  to  relate  zone  size  (Ax)  and  stress 
intensity  factor  (Kj).  On  the  other  hand 3 the  stress  field 
near  an  existing  c 20k  can  be  directly  related  to  the  crack 
propagation  mode t-  and  zone  size  used  in  the  numerical  solutions. 

The  stress  field  near  the  tip  of  a thin  crack  loaded  in 
tension  in  an  elastic  medium  (Reference  21)  is  given  by 


-7===-  Cos  (6/2) 
y2TTr 


1 - sin(0/2) sin(30/2) 
sin  (9/ 2)cos (30/ 2) 

1 + s in  (0/ 2)  sin  (30/ 2) 


where  the  following  sketch  defines  the  coordinate  system 


(6) 


It  is  important  to  note  that  Equation  6 is  based  on  the 
equations  of  elasticity  with  the  crack- free  surface  condition 
(no  tractions)  imposed  as  a boundary  condition.  The  numerical 
solutions  solve  exactly  the  same  equations  with  exactly  the 
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same  free  surface  boundary  conditions  on  the  crack  surfaces. 
Thus,  the  numerical  solutions  will  approximate  Equation  6.  The 
elastic  tensile  stress  causing  crack  propagation  along  the 
direction  of  the  existing  crack  (0=  0)  is  therefore  governed 
(near  the  crack  tip)  by 


a 


22 


K 


I 


( 2 tt  r ) 5 


(7) 


In  the  numerical  solutions,  crack  propagation  is  gov- 
erned by  the  magnitude  and  duration  of  tensile  stresses  in 
the  computational  cells  in  front  of  existing  cracks.  (In 
Section  4 , the  increase  in  tensile  stresses  in  front  of 
existing  cracks  in  the  numerical  solutions  was  demonstrated). 
If,  as  a first  approximation,  we  assume  that  the  typical 
distance  between  the  existing  crack  and  the  nearest  cell 
center  is  r = Ax/2,  and  if  we  assume  022  = for  crack  prop- 
agation in  the  numerical  solutions,  then  Equation  7 immediately 
relates  the  fundamental  stress  intensity  factor  to  the  para- 
meters oT  and  Ax,  i.e.. 


or 


(it Ax) ' 


(8) 


Note  that  Equation  8 completely  specified  the  effective 
stress  intensity  parameter  once  Oj  and  ax  are  specified. 

Thus,  each  zoning  and  variational  calculation  (with  icy  = 0) 
can  be  viewed  as  a variation  in  the  single  material  parameter 
Kj.  Table  5 indicates  the  effective  stress  intensity  factors 
for  the  three  cases  shown  in  Figure  42. 
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TABLE  5 


EFFECTIVE 

STRESS 

INTENSITY  FACTORS 
EQUATION  8 

ACCORDING 

Case 

* 

D 

Ax 

K** 

Ax 

(cm) 

* u 

psi  (cm) 2 

14 

10 

. 01 

1099 

1 

20 

. 005 

777 

15 

4 0 

.0025 

550 

All 

cases  use 

Ky-  0 , 

Oj  = 6500  psi,  v = 

.3 

**  Kj  for  ZnSe  from  Table  IV  is  600-1300  psi-inch'5 
which  is  960-2070  psi-cm'5 


Note  that  the  effective  Kj  used  in  the  numerical  calculations 
is  reasonably  close  to  the  experimental  values  for  ZnSe. 

An  important  verification  of  the  numerical  simulation 
of  Kj  based  on  Eq.  8 was  performed.  If  the  numerical  solu- 
tions do  approximately  simulate  a Kj  based  on  the  choice  of 
a.j.  and  Ax,  then  similar  crack  patterns  should  uevelop  for 
similar  values  of  Kj.  Figure  7,  page  15,  compares  Cases  15 
and  16  which  have  the  widest  variation  in  and  Ax  of  all 
the  cases  performed  in  this  study.  Table  6,  below,  summarizes 
the  relevant  parameters. 

TABLE  6 

EFFECTIVE  STRESS  INTENSITY  FACTORS  ACCORDING  TO 
EQUATION  8 


Case 

O 

(psi) 

D 

Ax 

Ax 

(cm) 

KI  , 

(psi-  (cm'5) 

15 

6500 

40 

.0025 

550 

16 

3000 

10 

.01 

507 
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Note  that  Oj  changes  by  a little  over  a factor  of  2 and  Ax 
varies  by  a factor  of  4 between  the  two  cases;  therefore, 

Kj  is  relat  ively  constant  ( ~ 8%  change). 

As  Figure  7,  page  15,  demonstrates,  the  crack  patterns 
for  the  two  cases  are  very  similar.  Both  targets  have  long 
cracks  extending  beyond  .05  cm  in  radius  (one  drop)  and  below 
a depth  of  .05  cm.  The  target  damage  in  the  two  cases  is  also 
very  close  (A/An  = 3.57  in  Case  15,  compared  to  3.36  in  Case 
16) . The  only  important  difference  in  the  two  patterns  is 
that  Case  16  (D/Ax  = 10)  has  developed  a crack  near  the  axis 
of  symmetry.  The  numerics  near  the  axis  of  symmetry  in  a rela- 
tively coarsely  zoned  problem  could  easily  cause  very  short 
tensile  oscillations  in  excess  of  the  low  a^.  value  of  3000  psi 
in  Case  16.  In  Case  15,  on  the  other  hand,  the  tensile  oscil- 
lations would  have  to  exceed  6500  psi  in  a very  finely  resolved 
numerical  solution.  For  this  reason,  we  would  discount  the 
cracks  near  the  axis  in  Case  15. 

There  is  still  much  to  be  understood  about  the  relation- 
ship between  the  material  physical  parameters  and  numerical 
solution  parameters  such  as  o^,,  k^,  and  zone  size.  However , 
the  re8ult8  to  date  indicate  that  the  numerical  crack  propa- 
gation model  used  in  this  study  can  be  systematically  related 
to  the  fundamental  stress  inten^kty  factor  Kj.  Additional 
efforts  should  be  directed  toward  improving  the  crack  initiation 
models  in  a form  which  is  independent  of  computational  cell 
size. 


< 
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